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Abstract 

Magnetic impurity adsorbed on one of the carbon planes of a bilayer graphene is studied. The 
formation of the many-body SU(2) and SU(4) resonances close to the bandgap is analyzed within 
the mean field Kotliar-Ruckenstein slave boson approach. Impact of enhanced hybridization and 

d 

'-^ . magnetic instability of bilayer doped near the Van Hove singularity on the screening of magnetic 

QJ ' moment is discussed. 
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I. INTRODUCTION 



Bilayer graphene (BLG) is a quasi-two-dimensional gapless chiral electron-hole system 
with parabolic bands with conduction and valence branches touching at a point The 
AB Bernal stacked bilayer has inversion symmetry. The non-zero bandgap can be induced 
by breaking this symmetry e.g. by applying transverse electric field. As opposed to conven- 
tional semiconductors, bandgap of BLG can be tuned, what is the reason of technological 
importance of this material |2j. The density of states (DOS) of biased BLG has the square- 
root singularities at the band edges, just like in one-dimensional systems. It is well known 
that electronic instabilities are expected at the crossing of the Fermi level (Ep) with Van 
Hove singularity (VHS). Thepossibility of magnetic ordering in bilayer graphene has been 
discussed in several papers 3|-[5j]. Since the chemical potential of BLG can be tuned one 
expects that formation of local magnetic moment can be also controlled by gate voltage. A 
strong enhancement of DOS at the Fermi level increases hybridization, what moves impurity 
state towards mixed valence (MV) regime or even suppresses the many body resonance. It is 
a challenging problem to discuss Kondo or mixed valence resonance under above mentioned 
conditions taking also into account possible magnetic instability in BLG matrix. Under- 
standing of the interplay of these competitive effects is important for prospect of Van Hove 
singularity engineering. 



II. MODEL 

Bilayer graphene [l| is composed of two coupled honeycomb lattices of carbon atoms. We 
consider AB-Bernal stacking, where the top layer has its A sublattice on top of sublattice B 
of the bottom layer. Indices 1 and 2 are used to label top and bottom layers, respectively. 
In the tight binding approximation the BLG Hamiltonian reads: 

U BLG = ( V2)4t + U 9 n%n%l - t £ (a^b aj(7 + h.c.) 

paicr paia 

-t ± J2( a ti^ + h.c), (1) 

i 

where a aia (b aio .) annihilates electron on sublattice A(B) in plane a = 1,2 at site i, (ij) 
represent a pair of in-plane nearest neighbors, p numerates sublattices, t, t± denote the in- 
plane and interlayer hopping integrals respectively, U g measures on-site graphene Coulomb 
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interaction and V g is the electrostatic bias applied to the two planes. The value of t (t = 2.7 
eV) is taken as the energy unit and we assume t± = 0.2, U g = 1 j^]. It is out of scope of 
the present paper to discuss the details of magnetic instabilities near VHS (see for e.g. 
||3]). For simplicity, following jlj] we adopt the mean field approach and define ferromagnetic 
ground state characterized by unequal magnetizations of the layers and different electron 
densities: 

where n and m denote electron density and magnetization per unit cell respectively and 
An, Am the corresponding differences between the layers. The analysis presented in this 
paper is addressed to TM adatoms located in the energetically favorable position -centre of 
a hexagon (hollow site). The symmetry of hollow site in BLG is Cs v and five d orbitals split 
into two orbital doublets corresponding to two-dimensional representation (E) and a singlet 
of one-dimesional representation A\. In Kondo physics the crucial role play electrons with 
unpaired spins. To get a microscopic insight into the orbital occupations of the concrete 



impurity ab initio calculations are required 



6|. In our model analysis based on Anderson 



Hamiltonian the picture is simplified by considering only the pair of orbitally degenerate 
states (S77(4)) or by discussing even simpler case when symmetry is further reduced to two 



spin-orbitals (577(2)) jg]. The impurity Hamiltonian for S77(4) symmetry reads: 

vW^C + k), (3) 



ma 



with \m\ = 1,2, e = -0.5. c ma = {1/VQ) Ei( eOTWa i«x + e tm ^ B b lla ) and (p iA{B ) denotes the 
angle between a fixed crystalline axes and the bond from site % to the impurity. To discuss 
the many-body problem we use mean field Kotliar-Ruckenstein slave boson approach [?| in 
the infinite impurity Coulomb interaction limit Uq — > 00. The finite Uq limit will be discussed 
elsewhere. Three introduced auxiliary boson fields e, p ma (er = ±1) project respectively onto 
the empty or single occupied state (with spin up or down). The slave boson Hamiltonian 
reads: 

7~Limp = ^ ^ ^ofmafma + V' ' ^ ^ { C ma ^ma fma + h.C.) 

ma ma 

+ A(7 - 1) + ^ma(Qma ~ f+Jma) (4) 
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The terms with Langrange multipliers A and A mcr are incorporated to prevent double oc- 
cupancy and to fulfill the complete relations. / = e + e + YlmaPmaPma, Qma = vturVma, 
Zrau = e + Pm*/(VQm*V~L ~ Qma) and pseudofermion operator f ma = z mu d ma . The slave 
boson parameters e, p a (p a = p ma = p_ TO(J ), A and X a are determined minimizing the ground 
state energy of (4). 



III. RESULTS AND DISCUSSION 



Fig. [U illustrates DOS of the mean field ferromagnetic state of BLG (g a (E)) for Fermi 
level located close to the upper gap edge. The spin splitting of near-gap conduction band 
states is visible. The low energy band structure presented in the inset shows a Mexican 
hat dispersion. The wave vector q measures the deviation from the centre of the valley 
K± = ±(4-7r/3a, 0), {k = K± + q). For the assumed positive potential applied to plane 1 
the wavefunction amplitudes of these states are higher for layer 1 than for layer 2. Even 
though the induced magnetic moments are extremely small, the spin polarization at the 
Fermi level Pblg — ^(e^ m ight be strong if E F is placed close to the upper edge. 
In the narrow energy range in this region the full spin polarization of electrons is observed. 
Apart from Van Hove singularities in the position of minima (maxima) of conduction band 
(valence band) also jumps of DOS are observed. They correspond to maximum (minimum) 
of energy in the centers of valleys (see the inset). Fig. |2] shows hybridization self energies 
for angular momentum symmetries m = 0,1,2. The presented hybridization function S m 
plays the role of self-energy of the impurity Green's function G (u) corresponding to (3), 
E m (u) = co — eo — Gq 1 (oj). As it is seen for hollow site hybridization of m = symmetry is 
strongly suppressed in the low energy range. The self-energies for m = 1 and m = 2 do not 
differ very much in this region, but for higher energies distinctive differences are observed 
and the dependences are strongly asymmetric with respect to the zero energy. Very roughly 
the electron part of m — 1 self energy corresponds to the hole part of m — 2. The observed 
splitting of the peaks around E — Ep pa 1.03 corresponding to M point of Brillouin zone are 
caused by interlayer coupling. The details of the spin dependent structure around band gap 
are presented on Fig. [HI where in addition to imaginary parts also real parts of self energies 
are plotted. Divergence of imaginary part for Ep located at the edge signals non-Fermi 
liquid behavior. The observed peaks and jumps of the real part substantially renormalize 
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FIG. 1: (Color online) Spin-resolved densities of states of BLG close to the bandgap plotted by solid 
(dotted) line for up (down) spin. Inset shows the corresponding Hartree-Fock bands, V g = 0.0185. 




FIG. 2: (Color online) Rescaled imaginary parts of hybridization functions for the impurity in the 
hollow site for three angular momentum symmetries T, m /V 2 (m = - dashed line, m = 1 - dotted 
line, m = 2 - solid line) compared with BLG density of states (dash dash dotted line). Inset shows 
the same hybridization functions as in the main picture, but in the narrower energy range. The 
assignment of the lines is preserved. 
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FIG. 3: (Color online) Real and imaginary parts of rescaled hybridization functions T,i a /V 2 plotted 
close to the bandgap. Solid and dotted lines represent imaginary parts for spin up and spin down 
respectively and dash dash dotted and dashed lines represent the real parts of the spin up and 
spin down self energies. For transparency imaginary part of hybridization function below the lower 
bandgap edge is extended four times. 

the effective impurity energy level in this region. We do not address here the subtle problem 



. We have 



of local moment formation in BLG, this has been recently discussed e.g. in 
only checked that in Hartree-Fock approximation the magnetic moment of impurities do not 
vanish for the assumed parameters. Figures H] and [5] show evolution of the expectation values 
of slave boson operators, magnetic moment and polarization of impurity versus Ep for 577(2) 
and S77(4) symmetries. The factors determining the character of the many-body resonance 
are the deepness of the atomic level with respect to Ep and hybridization strength. For the 
narrow energy range discussed, where the dramatic changes of the latter are observed, impact 
of the hybridization is dominant. The representative densities of states of the impurity 
(go tr (E)) are shown on Figures [61 Below the lower gap edge, where hybridization self-energy 
is small, structureless and spin independent (Figs. Ek), typical Kondo resonance is 
formed centered at Ep for SU(2) symmetry and shifted above Ep for S77(4) case (FigJB^). 
When Ep moves closer to the edge and both real and imaginary parts of self-energy are 
strongly enhanced broadening of many-body spectral function results for both symmetries 
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FIG. 4: (Color online) Evolution of SU(2) impurity state with the shift of the Fermi level in the 
narrow energy range around bandgap a) The expectation values of slave-boson operators e, p a and 
the impurity electron occupation n for coupling parameter V = 0.25. b) Polarization of impurity 
Po j impurity magnetization tuq and polarization of BLG (Pblg)- Insets of both pictures present 
the corresponding dependences for V = 0.1 for Ep located above the upper gap edge. 




FIG. 5: (Color online) The same as in Fig. 4 but for SU(4) symmetry. 
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FIG. 6: (Color online) Representative impurity densities of states a) V = 0.25, E F located below 
lower bandgap edge (unpolarized BLG): SU{2) E F rj 0.48580 dashed line, SU(2) E F ps 0.49080 
dash dash dotted line, SU(4) E F « 0.48580 dotted line, 5t/(4) £ F ps 0.49080 solid line, b) 

V = 0.25, DOS of orbital MV state SU(A) -»• 5i7(2) for £ F « 0.50909 (fully spin polarized BLG), 
c) V = 0.25, SU(2), Ep located between the upper bandgap edge and nearby step in BLG DOS 
(polarized BLG) (E F ~ 0.50913), solid line for spin up, dotted line for spin down, d) the same 
as (c), but for SU(4) symmetry, e) the same as (c), but for V = 0.1, f) the same as (d), but for 

V = 0.1. 

and for 577(4) additional structure at the edge is also visible reflecting a new pole in the 
impurity Green's function occurring close to the edge, but still in the region of finite density 
of states of BLG. Within the gap a resonance can be also formed, and in this case the 
resonance is not masked by finite DOS and takes the form of Dirac delta like line (Figs. 
Eb,d). Of interest are also dips in the impurity DOS (Figures Eb-e), occurring for energies 



S 



where steps in the imaginary parts and peaks in real parts of self energies are observed (Fig. 
[3]). The narrow energy region above the upper edge is of special interest for spintronics, 
because spin polarization of BLG matrix results in spin splitting of resonances (Figs. Efc-f). 
For typical coupling parameter corresponding to the equilibrium distance of impurity from 
the graphene layer (assumed value V = 0.25 corresponds to Co impurity-graphene coupling 
js]), the calculated average occupation considerably deviates from n = 1 i.e. impurity is 
in mixed valence state for both symmetries (Figs. |6b,d, Figs. H^JS^i). More interesting is 
the case of a weaker coupling, when the spin polarized Kondo state is formed (Figs. [6fe,f). 
Possible reduction of the coupling can be achieved e.g. by chemical engineering (ligands), 
or introducing ultrathin dielectric film separating impurity from graphene layer, or putting 
the adatom further from the surface by STM methods |IJ] . For Fermi level located between 
the upper edge and nearest higher in energy step in BLG density of states the impurity 
polarization for SU(2) symmetry is opposite to BLG polarization, whereas for SU(4) has 
the same sign (Figs. HbJSb)- This tendency is observed both in MV and in Kondo states 
and is in agreement with earlier calculations for regular density of states of the electrodes 
performed with the use of different techniques 10h12|]. This behavior reflects the fact that 
despite the same sign of the spin splitting for both symmetries (opposite to BLG spin 
splitting), the dominance of spin resolved densities of states at the Fermi level is different 
due to significantly different location of unperturbed many-body peaks in 577(2) and 577(4) 
cases. Interesting observation is a possibility of full spin-polarization of impurity state in 
the region of full polarization of BLG matrix and reverse of magnetic moment when Ep 
enters the region of full BLG spin polarization. Since the Fermi level can be shifted by gate 
voltage this gives the way of electric control of magnetic moment, what is a highly desirable 
property for spintronics. 

In summary, present investigation only signals the interplay of different factors on the 
many-body physics of magnetic impurity in bilayer graphene doped near the Van Hove 
singularity at the bandgap edges. The exact description of the crossing of the Fermi energy 
with a Van Hove singularity in the DOS would require more sophisticated methods, both in 
defining the magnetic state of the bilayer going beyond mean field scheme and in discussion 
of screening processes, where analysis outside SBMFA is required. In addition to hole- 
like processes discussed above also electron-like charge fluctuations might play the role in 
reconstruction of Kondo resonance and this topic will be the subject of a forthcoming paper. 
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